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ABSTRACT 

We investigate the relationship of the H II region luminosity function (H II LF) to the H II region size 
distribution and density wave triggering in grand-design spiral galaxies. We suggest that the differential 
nebular size distribution is described by a power law of slope ~ —4, with flattening at radii below ~ 130 
pc. This contrasts with the conventional exponential description, but it is physically and quantitatively 
^> ' consistent with the typical observed value of —2 for the HII LF slope. 

\ We have developed an interactive code that computes spatial isochrones for the evolving loci of spiral 

CO i density waves in disk galaxies. This allows comparison of the nebular spatial distribution with the spatial 

isochrones for simple rotation curve parameters. Our comparisons for four grand-design galaxies suggest 
, that the corotation radius Tco coincides with the outer ends of the star- forming arms. This value for 

Tco yields the best spatial correspondence between the H II regions and the isochrones, and also appears 
to yield a coincidence between the Inner Lindblad Resonance with the radial onset of star formation in 
the arms. Thus, we suggest that isochrones offer a new, simple, and effective technique for determining 
' Tco, and thus the spiral pattern speed. However, application of the isochrones also demonstrates that 

evolution of the nebular population is difficult to spatially isolate in these galaxies. 

O ' Subject headings: H II regions — galaxies: fundamental parameters — galaxies: kinematics and 

^ , dynamics — galaxies: spiral — ISM: kinematics and dynamics 



O 

cn 
o 



>< 



1. INTRODUCTION tion in A^*, the number of ionizing stars per object: 

The H II region luminosity function (HII LF) is an im- N{N ) dN ~ B dN (2) 

lortant diagnostic of star formation properties in galax- * 

M ; ies. It is a fairly straightforward parameter to deter- We then performed Monte Carlo simulations, in which we 

5^ , mine for nearby galaxies (e.g., Kennicutt, Edgar, & Hodge draw from this distribution, and also draw the indi- 

1989; Banfi etal. 1993; Rozas, Beckman, & Knapen 1996a; vidual stellar masses from a Salpeter (1955) initial mass 

Thilker etal. 2002), and is usually presented in the differ- function (IMF). These models showed that reported vari- 

ential form as a power law: ations of the measured slope a (equation 1) in different 

environments do not require varying (5, the slope of the 

NiL) dL — A L^" dL (I) ^* distribution (equation 2), but rather, can be explained 

' by variations in the star formation history and the up- 
per cutoff in the N^, distribution. An important feature is 

where N{L) dL is the number of nebulae with Ha lumi- that stochastic effects for objects with small N^, cause a 

nosities in the range L to L -I- dL. Occasionally, the cu- flattening in the HII LF below a certain luminosity igat- 

mulative H II LF is given, but we note that the differential Sparse clusters subject to this effect are termed "unsatu- 

form is a more powerful diagnostic of the underlying prop- rated" with respect to the IMF, and clusters with higher 

erties (Cioffi & ShuU 1991). V, we call "saturated". We identified igat with the Ha 

In an earlier paper (Oey & Clarke 1998, hereafter Pa- luminosity associated with the most massive star in the 

per I), we modeled the behavior of the HII LF, and iden- IVIF, although we refer readers to Paper I for a complete 

tified ways in which it can be used to examine the current discussion. 

properties of massive star formation and recent (^ 10 Myr) A basic test on the nature of the HII LF is its relation 

star formation history. We assumed a power-law distribu- to the size distribution of the H II regions. Although the 

^Participant in the 1999 STScI Summer Student Program. 

^Present address: University of Colorado, Dept. of Aerospace, Engineering Sciences, Boulder, CO 80309 
■^Present address: University of Florida, Dept. of Astronomy, P.O. Box 112055, Gainesville, FL 32611 
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H II LF has empirically been described by a power law, the 
corresponding size distribution has been described by an 
exponential law (van den Bergh 1981; Hodge 1983, 1987). 
As discussed above, Paper I found that the slope breaks 
and other features in the HII LF are nevertheless best 
explained by a simple power law in A^* (equation 2 of Pa- 
per I). Therefore we should naively expect the size distri- 
bution to also follow a power-law distribution. In §2 below, 
we examine whether observations of the H II region size 
distribution are consistent with a simple power-law rela- 
tionship. 

Another test on the nature of the H II LF addresses its 
evolutionary properties. The models constructed in Pa- 
per I included the evolution of the H II LF resulting from 
a single burst creation scenario. In Figure 1, we show rep- 
resentative zero-age and 7 Myr models (Figures 2a and 
c of Paper I). The zero-age model shows the flattening 
below isat- As time goes by, the entire HII LF shifts 
to lower L. However, when considering a fixed range of 
L, this can result in an apparent elimination of the slope 
break in strongly evolved populations, due to the differ- 
ential fading effect. If this model turns out to be correct, 
it would then be possible to distinguish, for example, be- 
tween an older nebular population where high-luminosity 
objects were created in a single burst, and continuous cre- 
ation of low-luminosity H II regions. 

Paper I applied these results to observations of arm and 
interarm regions of grand design spiral galaxies, suggesting 
that arm populations represent the current burst of star 
formation, while interarm objects represent an aged pop- 
ulation remaining in the wake of the spiral density waves. 
Observations in the literature sometimes suggest that in- 
terarm populations show a slightly steeper H II LF than 
the arm populations (e.g., Kennicutt etal. 1989; Rand 
1992). This phenomenon can be explained by the two- 
slope structure in the zero-age models, which will cause a 
slope measurement to be flatter when fltted over L below 
Lsai (see Figure 1). In contrast, the evolved populations 
show the intrinsic slope of the high-L regime. Figure 5 
of Paper I shows that this interpretation of evolution be- 
tween the arm and interarm H II regions is consistent, in 
all cases, with available observations of six grand design 
spiral galaxies. 

Observations of the H II LF in the literature to date have 
been subjectively and coarsely binned into the arm and in- 
terarm categories. A more powerful test of this evolution- 
ary scenario would be to rebin the H II regions more quan- 
titatively into regions corresponding to fixed age ranges be- 
hind the leading edge of the spiral arms, as inferred from 
the spiral arm pattern speeds. We might then be able to 
examine the competing contributions in interarm regions, 
of aging nebular populations versus bona fide interarm star 
formation. These are our motivations for carrying out the 
studies described in §3. As an important by-product of 
this study, we developed a technique that compares nebu- 
lar positions with spatial isochrones for the spiral density 
waves. This offers a new method for locating the corota- 
tion resonances in spiral galaxies, and is described in §4. 

2. THE FORM OF THE H II REGION SIZE DISTRIBUTION 

As mentioned in the Introduction, the conventional pa- 
rameterization of the H II region size distribution as an 
exponential (e.g., van den Bergh 1981; Hodge 1983, 1987) 



is intuitively at odds with the power-law form of the HII 
LF. Since the HII LF's behavior appears to be fairly well- 
understood (e.g.. Paper I), it is important to reexamine the 
size distribution within the same framework, since they are 
physically related. In particular, we might also expect the 
size distribution to follow a power law of the form, 

Nr{R) dR DC R-'' dR , (3) 

where R is the nebular radius. This has indeed been sug- 
gested for some galaxies (e.g., Pleuss, Heller, & Fricke 
2000). 

Since the nebular Ha luminosity L is integrated from 
the nebular volume emission, we should expect L (x R? 
and ^ a i?^. Therefore, the differential size distribution 
should be: 

Nr{R) dR = N{L) ^ dR(x R'^'^" dR . (4) 
dR 

Hence we find that the slopes a and b of the H II LF and 
size distribution, respectively, are related as 6 = 2— 3o. For 
the usual value a = 2 found for the high-luminosity, "sat- 
urated" H II regions, we therefore predict a corresponding 
size distribution having 6 = 4. 

As mentioned earlier, a turnover in the H II LF caused by 
stochastic effects for "unsaturated" objects occurs at a lu- 
minosity Lgat • Since igat is associated with the most mas- 
sive stars in the IMF, their Stromgren radius i?sat should 
be a corresponding inflection point in the size distribution. 
This flattening at small R can mimic an exponential dis- 
tribution. For the most massive stars, stellar atmosphere 
models suggest logLgat/ ergs"-'^ ^ 38.0 to 38.5 (Schaerer 
& de Koter 1997; Vacca, Garmany, & ShuU 1996; Panagia 
1973), implying i?sat ~ 130 pc for a density of 1 cm~^. 
This estimate for iJgat should be taken as an upper limit, 
since the density is likely to be higher in the actual, central 
star forming area. 

We can now compare observations of the H II region 
size distribution to our predictions for slope 6 = 4, with a 
flattening in R below logi?sat/pc ^2.1. We refltted the 
nebular diameter distributions with power laws for a num- 
ber of galaxies in the literature. Since size determinations 
can be sensitive to different methods, and thus to ensure 
relatively consistent comparisons, we used data from only 
two groups of investigators: the lAC Hertfordshire collab- 
oration (Knapen et al. 1993; Rozas et al. 1996b; Knapen 
1998), and Youngblood & Hunter (1999). Hodge's group 
(e.g., Hodge etal. 1994); Miller & Hodge 1994; Hodge 
etal. 1999) and Petit's group (e.g., Petit 1993, 1996, 1998) 
have contributed many studies of H II region size distribu- 
tions; however, their target galaxies generally do not of- 
fer enough objects to adequately sample the nebular size 
distribution at logD/pc > 2.4, where D is the nebular di- 
ameter. Likewise, we include only the three galax:ies from 
the Youngblood & Hunter (1999) sample that provide ad- 
equate sampling in this regime. 

Figure 2 shows the size distributions for the nine galax- 
ies listed in Table 1. We fitted power law slopes for 
\ogN{D)/pc > 2.3, with the bins weighted by ^/W. The 
lAC-Hertfordshire galaxies are grand design spirals and 
the Youngblood & Hunter galaxies are Magellanic irregu- 
lars. Both these sets of galaxies have H II LFs extending to 
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Fig. 1. — Model HII LFs: (a) zero-age; and (6) evolved single-burst, from Oey & Clarke (1998; Figures 2a and c). The dotted line marks 



high luminosities, consistent with no upper cutoff in lumi- 
nosity. Our fits to the size distributions therefore should 
not be influenced by possible truncation at the upper end 
of the HII LF. It appears that the H II region diameter 
distributions shown in Figure 2 are consistent with power 
laws in the given regime. Objects with smaller diameters 
show a flattening in the size distribution, qualitatively con- 
sistent with the flattening in the HII LF for statistically 
unsaturated objects. 

The small dynamic range and small number statistics 
cause some ambiguity and large uncertainties in the fitted 
slope — fe, as shown in Figure 2. We also caution that we 
cannot rule out the existence of systematic effects. Nev- 
ertheless, the average slope, weighted inversely by the un- 
certainties, is (6) = 4.12 with a standard deviation of 0.90. 
This preliminary result agrees with the value predicted 
above of 6 = 4. Figure 2 also shows that the turnover in 
the diameter distribution is also consistent with the pre- 
dicted value of logi^/pc ~ 2.4. We therefore suggest that 
the H II region size distribution is intrinsically described by 
a power-law function in the regime beyond the turnover, 
consistent with its physical relation to the HII LF. While 
consistent functional forms are expected for the size distri- 
bution and HII LF, further confirmation of the empirical 
power law for the former is desireable, with attention to 
systematic uncertainties in size determinations. 

3. TESTING FOR AN EVOLVING HII LF 

Based on Paper I, we approximate that the H II regions 
found along the leading edge of spiral arms represent the 
youngest, zero-age population of objects. H II regions be- 
hind the leading edge therefore represent a somewhat older 
population remaining in the wake of the spiral density 
wave. If the galactic rotation curve and spiral arm pattern 
speed are known, it is in principle possible to estimate the 
ages of objects based on their location with respect to the 



leading edges of the spiral arms. This would also assume 
that the velocities of the H II regions do not deviate signif- 
icantly from the galactic rotation curve during the passage 
of the density wave. 

We have developed a code that computes spatial 
isochrones from the rotation curve, given the current lo- 
cus of the spiral arms. We defer discussion of the code and 
spiral arm kinematics until §4 below, and here use the spa- 
tial isochrones to search for evolutionary features as pre- 
dicted in Paper L We examined four of the grand design 
spiral galaxies studied in Paper I, while the two remain- 
ing galaxies included in that paper, NGC 3631 and NGC 
6814, did not have a simple two-arm structure suitable 
for this analysis. Figure 7 below shows Ha images of the 
galaxies: M51 (Rand 1992), MlOO (Knapen etal. 1993), 
NGC 157, and NGC 6951 (Rozas etal. 1996). Five sets of 
spatial isochrones are calculated for each of the two dom- 
inant arms in these galaxies, calculated at intervals of 12 
Myr. The color scheme for times -12, 0, 12, 24, and 36 
Myr are: red, green, yellow, yellow, red for one arm; and 
blue, yellow, green, green, blue for the other. The zero- 
age isochrones correspond to an actual spatial fit to the 
arm positions. The adopted corotation radius is shown 
with the black dashed line, and the inner and outer Lind- 
blad resonances are shown in the solid and dotted black 
lines, respectively, assuming m = 2 spiral arms for these 
galaxies. 

Unfortunately, it appears that the isochrone spatial 

separation that is sensitive to the nebular evolutionary 
timescale is of order the spatial scatter of the objects along 
the spiral density wave. Figures 7a — c show that the 
best-defined arms are enclosed within the -12 to +12 Myr 
isochrones, but show a spatial scatter that precludes fur- 
ther distinction. NGC 6951 is an exception (Figure 7d), 
showing larger separation between the isochrones that also 
corresponds to the nebular spatial distribution. 
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Table 1 

References for galaxy diameter distributions 



Galaxy 

MlOO 
NGC 157 
NGC 1156 
NGC 2366 
NGC 3631 
NGC 4214 
NGC 6814 
NGC 6764 
NGC 6951 



Type- 



Reference 



-II 



Sc(s)I 
Sc{s)I- 
SmlV 
SBmlV-V 
Sbc(s)II 
SBmlll 
Sbc{rs)I-II 
SBbc(s)'' 
Sb/SBb(rs)I.3 



Knapcn (1998) 



Rozas et al 
Youngblood 
Youngblood 
Rozas et al. 
Youngblood 
Knapen etal. (1993) 
Rozas et al. (1996b) 
Rozas et al. (1996b) 



(1996b) 

& Hunter (1999) 
& Hunter (1999) 
(1996b) 

Hunter (1999) 



-From Sandagc & Tamann 1987. 
''Type as listed by Rozas et al. (1996b). 
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Differential H II region diameter distributions for galaxies in Table 1, fitted with power laws for logD/pc > 2.3. The fitted 
and standard deviation are shown [note that the fitted slope for logAf(logD) vs. logD is 1 — b ]. Empty bins are shown as 
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Figures 3-6 show the H II LF for sub-groups of H II re- 
gions binned within isochroncs for 12 to +12 Myr (panels 
a); to 24 Myr (panels b), and 12 to 36 Myr (panels c). As 
discussed above from Paper I, the H II LF naturally shows 
a break in slope below the high-luminosity power-law tail, 
that is caused primarily by the transition to unsaturated 
stellar clusters dominated by small-number statistics. For 
older nebular populations, this slope transition at Lgat is 
expected to shift toward lower luminosities, while the high- 
luminosity power-law slope remains the same (Figure 1 
above; Paper I). The vertical dotted line in these figures 
shows the luminosity bin whore wo estimate Lgat to occur 
for the youngest binning; this value is plotted for reference, 
and remains constant between the three HII LF subsets 
for each galaxy. 

It is apparent in Figures 3-6 that the location of Lgat 
is essentially indistinguishable between the three subsets. 
The transition point generally appears to be more am- 
biguous in the oldest binnings (panels c) as compared to 
the youngest (panels a), but more galaxies must be exam- 
ined to determine whether that pattern is real. If so, then 
it would be consistent with a greater contribution from 
evolved populations (Paper I). The ambiguities again re- 
flect our discussion above, that the isochroncs for these 
galaxies are insensitive to aging effects, since the spatial 
scatter is of order the separation for isochroncs over the 
nebular evolution time scale. 

A number of other factors probably also contribute to 
the insensitivity of the isochroncs to aging efi'ects. First, 
the models in Paper I (e.g.. Figure 1 above) assumed 
the instantaneous creation of all objects in the HII LF, 
whereas triggering by the spiral density wave most likely 
has a finite timescale. Since most H II regions recombinc 
on timescales of ^ 10 Myr, we would not expect any sig- 
nificant population to be visible at 24 Myr. For the most 
part, this is consistent with the fairly negligible Ha emis- 
sion at the 24 Myr isochrone for most of the arms in M51, 
MlOO, and NGC 157. NGC 6951 shows spatially more ex- 
tensive star formation, with H II regions also seen along 
the 24 Myr isochrone; this therefore implies an extended 
timescale for star formation at any given location. An 
extended star-formation timescale is also equivalent to a 
spatial scatter for objects of any given age, as mentioned 
above. For the HII LFs, this will blur the distinction be- 
tween the youngest subset and next-to-youngest. It may 
also be that the luminosity fading of the H II regions is 
even faster than our extreme model in Paper 1 (see Pa- 
per I, Figure 2/). Another factor that affects the quan- 
titative comparison to the H II LF evolution is the uncer- 
tainty in identifying and fitting the locus of the zero-age 
density wave. Finally, our simple analytic description for 
the rotation curve and pattern speed are idealized forms 
that do not describe the details of the gas kinematics (see 
§4 below), which also affect the location of the spatial 
isochroncs. 

Thus, our test of the evolutionary features in the HII 
LF for these galaxies is inconclusive, though an upper limit 
on the evolutionary effects is indicated. For most of the 
galaxies, the spatial isochroncs do not offer the necessary 
separation to isolate populations of H II regions triggered 
in the arm density waves. 



4. LOCATING COROTATION WITH SPATIAL ISOCHRONES 

We now describe our analysis of spatial isochroncs, 
which offer some interesting insights on spiral arm kine- 
matics. 

While measurements of rotation curves are fairly 
straightforward, historically it has been much more dif- 
ficult to estimate the spiral density wave pattern speed, or 
equivalently, to determine the location of the corotation 
radius. Several methods for determining the corotation ra- 
dius or pattern speed have been proposed in the past. The 
earliest proposal was based on Lin's (1970) conjecture that 
organized, two-armed spiral patterns exist only within the 
corotation circle, and that arms become more fragmented 
outside. Shu, Stachnik, & Yost (1971) were among the 
first to apply this conjecture to physical galaxies, placing 
the corotation where the distribution of H II regions is 
seen to end. Tremainc & Weinberg (1984) developed an 
approach for pattern speed determination based on mea- 
surement of surface brightness and velocity distributions 
of luminosity density tracers, assuming continuity in the 
tracer. Canzian's (1993) method for determining corota- 
tion relies on the qualitative difference in the residual ve- 
locity field inside versus outside of corotation. Elmegreen, 
Elmegreen, & Seiden (1989) identified resonances through 
breaks in spiral arms, and Cepa & Beckman (1990) utilized 
the radial star formation efficiency. Elmegreen, Wilcots & 
Pisano (1998) determined corotation from the transition 
of the inward to outward streaming gas velocity. Finally, 
Clemens & Alexander (2001) used high-resolution spectro- 
scopic data to find signatures of the pre- and post-shock 
gas across the spiral arms; this information determines the 
shock velocity, and hence, the pattern speed. In what fol- 
lows, we explore a new approach that uses the nebular 
spatial distribution to constrain the corotation radius. 

For a galaxy with a flat rotation curve, the corotation 
radius Tco is related to the angular pattern speed ftp as, 



where Wcirc is the circular velocity. Assuming that star 
formation is triggered by spiral density waves, the locus of 
the youngest H II regions essentially indicates the location 
of the density wave. As the nebulae age, their positions 
will be determined by both the spiral arm pattern speed 
and the galaxy's rotation curve. It is therefore possible to 
determine spatial isochroncs for a given galaxy, that repre- 
sent the positions of objects that all formed simultaneously 
within a density wave. 

Our code computes spatial isochroncs, given an initial 
locus for the spiral arm pattern. The H II region positions 
are deprojected and converted to polar coordinates, and 
we fit the spiral pattern at the arm's leading edge, using a 
polynomial function. We then compute spatial isochroncs 
from the locus of the spiral arm pattern, rotation curve, 
and putative corotation radius. We assume a flat rotation 
curve, except for the innermost regions r < rdic, where we 
assume v ocr (Table 2). 

Our code includes a widget to easily adjust the assumed 
corotation radius. Varying rco by hand, we found a value 
of Tco that determines a set of isochroncs having a notice- 
able correspondence to the spatial configuration of the H II 
regions. These values of Tco turn out to correspond to the 
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Fig. 3. — 6. — HII LFs for objects binned between isochrones for M51, MlOO, NGC 157, and NGC 6951. Panels a,b and c respectively 
show Bins 1-2 (-12 to +12 Myr), Bins 2 - 3 (0 - 24 Myr), and Bins 3 - 4 (12 - 36 Myr). 
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ends of the spiral arms in the star-forming disk. As we 
discussed above, other investigators (e.g., Slni et al. 1971; 
Rohlfs 1977; Kenney etal. 1992; Zhang, Wright, Alexan- 
der 1993) have identified r^o with the ends of the spiral 
arms in the past. For three of our four galaxies. Figure 7 
shows that these values of Tco also yield inner Lindblad 
resonances (ILRs) that coincide well with the radial on- 
set of significant star formation in the arms. Figure 7b 
also suggests a correspondence between the outer Lind- 
blad resonance (OLR) and the edge of the star-forming 
disk in MIOO. 

For all four galaxies, the regions near the ILRs dis- 
play the characteristic "twin-peaks" morphology (Kenney 
etal. 1992) or "twin-peaks-plus-inner-bar" (Zhang et al. 
1993), together with an annular zone between the twin 
peaks and the inner bar /nucleus that is swept clear of star 
formation. This combination of central morphologies is 
common among mid- to early-type galaxies, and is related 
to the orbital behavior of gas near the double inner Lind- 
blad resonances that are present in these galaxies with 
high central potential wells (Zhang et al. 1993, and refer- 
ences therein) , namely, the inner- inner-Lindblad resonance 
(IILR) and the outer- inner-Lindblad resonance (OILR). 
This morphology allows the effective central channeling of 
gas and fueling of nuclear activity. Note that our sim- 
plified assumed rotation curves yield only one ILR which 
presumably is located between the IILR and OILR of the 
actual rotation curves. 

For M51 and MIOO, Figure 7 shows that setting Vco 
to the tips of the spiral arms causes the isochrones to 
bunch or converge, but not cross over, at azimuthal lo- 
cations near large clumps of nebulae. These clumps are 
not confined within the isochrones, and extend to radii 
substantially beyond them; however, their azimuthal coin- 
cidence with the isochrone convergences strongly suggests 
that these loci are associated with dynamical triggering of 
the star formation. These patterns are most likely related 
to non-c;ircular motions and asymmetry in the mass distri- 
butions. Indeed, both M51 (see Garcia-Burillo et al. 1993) 
and MIOO (Garcia-Burillo et al. 1998) show evidence of an 
inner tri-axial mass distribution. It is also apparent in 
Figure 7 that the spatial isochrones only cover a fraction 
of the active star-forming regions in the disks, which again 
must be caused by non-circular kinematics. However, our 
experimentation with rotation curve parameters showed 
that to first order, our simple assumptions for the rota- 
tion curve (flat, with a linear inner portion) worked best 
in matching the Ha morphology of the arms. 

A counterexample is shown in Figure 8 for M51. 
One previous estimate for Tco in this galaxy is 9.6 kpc 
(Elmegreen etal. 1992). It is apparent that the isochrones 
do not correspond to the H II region positions as cleanly as 
in Figure 7a. The arms show a smooth, continuous mor- 
phology in the region around the suggested rco, whereas 
the predicted isochrones change character and diverge. 

Table 2 shows our assumed rotation curve parameters 
and our estimates for the rco, ILR, and OLR. We now 
compare our results for rco with those in the literature, 
in particular, estimates derived from kinematic models. 
There are no modeled estimates of rco for NGC 6951. 



4.0.1. M51 

Our analysis suggests that the corotation radius is twice 
as large as most previous determinations for M51. Both 
morphological estimates (e.g., Vogel et al. 1993; Elmegreen 
etal. 1992; Knapen et al. 1992) and kinematic model es- 
timates (Garcia-Burrillo et al. 1993b) found rco of 126" — 
161", in contrast to our estimate of 301". We note that 
M51 is one of the few examples where the arm definition 
is narrow enough that the i?-band and Ha or CO observa- 
tions can clearly distinguish the stellar arms from the new 
star-formation (Rand & Kulkarni 1993; Garcfa-Burrillo 
etal. 1993a). At corotation, the youngest regions should 
switch from the inside edge of the arm to the outside; yet 
the data shown by Rand & Kulkarni (1990) and Garcfa- 
Burrillo (et al. 1993a) show no evidence for this transition 
at the smaller rco suggested in the past. 

4.0.2. MIOO 

For MIOO, our corotation radius is larger than previous 
estimates, but by a smaller margin. Estimates by Garcfa- 
Burrillo et al. (1998) and Canzian & Allen (1997), based on 
kinematic modeling, and also Elmegreen etal. (1992) all 
agree fairly well with the range of 82" — 114" estimated for 
rco by Sempere et al. (1995). Our estimate of Vco ~ 154" is 
larger, but is based on more simplistic morphological and 
kinematic assumptions; most of the above studies consider 
the more realistic, bar-dominated nature of the rotation 
curve. We also note that the ends of the arms are not 
clearly delineated in this galaxy. 

4.0.3. NGC 157 

Our estimate of 

^co '"^ 56" agrees well with estimates 
in the literature based on detailed kinematic models. 
Sempere & Rozas (1997) find rco ~ 50", and Fridman 
etal. (2001) find 36" < rco < 63" by two different meth- 
ods. Elmegreen et al. (1992) also estimate corotation to be 
at 56" based on morphological arguments. 

Spatial isochrones thus appear to offer a new technique 
for constraining the corotation radius for a first-order ap- 
proximation. We can set a spatial correspondence be- 
tween: (1) Corotation and the end of the nebular spiral 
arms; (2) Isochrone bunching and H II region clumps; and 
possibly (3) Outer Lindblad resonance and the edge of 
the star-forming disk. These physical coincidences sug- 
gest that the location of corotation has been identified to 
first order. Thus, while some authors often identified the 
ends of the star-forming arms with the corotation reso- 
nance, and others with the OLR, our analysis supports 
the former. We also note that in wide-format, blue photo- 
graphic images available via the NASA/IPAC Extragalac- 
tic Database (e.g., from the Sandage-Bedke Atlas), there 
is no evidence in any of our galaxies that the dust lane 
crosses over from the inside to the outside of the arms at 
any point inside our adopted rco- 

5. CONCLUSION 

To better understand the behavior of the H II LF and its 
ability to probe global star-formation properties, we have 
re-examined its relation to the H II region size distribu- 
tion, and investigated its evolution with respect to density 
wave triggering in spiral galaxies. 
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Fig. 7. — H« images of the sample galaxies, overlaid with spatial isochroncs for —12, 0, 12, 24, and 36 Myr. These are respectively shown 
as red, green, yellow, yellow, red for one arm; and blue, yellow, green, green, blue for the other. The black dashed line shows the adopted 
corotation radius; the inner and outer Lindblad resonances are shown by the solid and dotted lines, respectively. Figure 7(a) shows M51, 
with Ha data from Rand (1992). 
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Fig. 7.— 



(6) MlOO; Ha data from Knapon et al. 1993 
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Fig. 7.— 



(d) NGC 6951; Ho data from Rozas et al. 1996 
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Fig. 8. — Same as Figure 7a, but with the M51 corotation at 9.6 kpc suggested by Elmegreen etal. 1992. 
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Table 2 

Galaxy and resonance parameters 



Galaxy 


Distance 




^circ 


^CO 


ILR 


OLR 


Reference'' 




(Mpc) 


(kms-l) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 




M51 


9.6 


220 


0.93 


14 


4.1 


24.9 


Rand (1993) 


MlOO 


16.1 


210 


1.56 


12 


3.5 


20.5 


Scmperc etal. (1995) 


NGC L57 


22.2 


170 


0.51 


6 


1.8 


10.2 


Sompcrc &■ Rozas (1997); Rvdor ctal. (1998) 


NGC ()l)51 


19.0 


2:i") 


0. 10 


10 


2,1) 


17,1 


«uzas el al, (2002) 



^References for rotation curve parameters; rdrc is the radius beyond which the circular velocity is taken to be constant, 
at value Ucirc- 



We have analytically re-evaluated the relationship be- 
tween the HII LF and the H II region size distribution. 
We suggest that the size distribution is intrinsically de- 
scribed as a power law, as physically related to the HII 
LF, rather than its conventional description by an expo- 
nential relation. We find that the physical relation be- 
tween nebular luminosity and size implies that the slopes 
of the size distribution b and the luminosity function a 
should be related as & = 2 — 3o,. The normally observed 
value a = 2 thus yields 6 = 4, which we find to be consis- 
tent with observations. We find that the slope transition 
radius logi?sat ^ 2.1 pc corresponding to the transition 
luminosity logigat ^ 38.5 is also consistent with observa- 
tions. 

We have developed a code to compute spatial isochrones 
from the rotation curve, given the current locus of the spi- 
ral density waves. Our code allows quantitative compar- 
ison of the nebular spatial distribution with the spatial 
isochrones. We find that the best correspondence for the 
arm H II regions occurs when the corotation radius coin- 
cides with the ends of the star-forming arms. This value 
for Tco also yields a coincidence between the Inner Lind- 
blad Resonance with the radial onset of star formation in 
the arms for three of the four galaxies. There may also 
be a hint that the Outer Lindblad Resonance coincides 
with the outer edge of the star-forming disk in at least one 
galaxy. The spatial isochrones therefore appear to offer a 
new technique for determining the corotation radius and 
principal resonances. 

We had hoped that these isochrones wotdd isolate sub- 
sets of H II regions that could then be used to study 
evolutionary effects of the nebular population. However, 
for three of the four galaxies, the isochrones do not yield 
adequate spatial separation over the nebular evolution 
timcscalc, relative to the spatial scatter at any given age. 
Encouragingly, we do find that the spatial concentration 
of the objects along the isochrones is consistent with the 
overall nebular evolutionary timescale in these galaxies. 
However, the data do suggest periods of extended star for- 
mation at any given location that blur the diagnostics. 
Combined with uncertainties in the kinematic parameters, 
this greatly limits the ability to obtain clear quantitative 
constraints on H II LF evolutionary parameters as related 
to the density wave passage. The observed rate of evolu- 
tion, however, does set limits on the importance of this 
process in producing the observed arm-interarm variation 
in the HII LF. 
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